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ABSTRACT 

Spectral lag, which is defined as the difference in time of arrival of high and low energy pho- 
tons, is a common feature in Gamma-ray Bursts (GRBs). Previous investigations have shown 
a correlation between this lag and the isotropic peak luminosity for long duration bursts. How- 
ever, most of the previous investigations used lags extracted in the observer- frame only. In this 
work (based on a sample of 43 Swift long GRBs with known redshifts), we present an anal- 
ysis of the lag-luminosity relation in the GRB source-frame. Our analysis indicates a higher 
degree of correlation —0.82 ± 0.05 (chance probability of ~ 5.5 x 10 -5 ) between the spectral 
lag and the isotropic peak luminosity, L- lso , with a best-fit power- law index of —1.2 ± 0.2, 
such that L iso oc lag -1,2 . In addition, there is an anti-correlation between the source-frame 
spectral lag and the source-frame peak energy of the burst spectrum, £' p k(l + z). 

Key words: Gamma-ray bursts 


1 INTRODUCTION 

Gamma-ray Bursts (GRBs) are extremely energetic events and pro- 
duce highly diverse light curves. A number of empirical correla- 
tions between various properties of the light curves and GRB en- 
ergetics have been discovered. However, the underlying physics of 
these correlations is far from being understood. 

One such correlation is the relation between isotropic peak lu- 
minosity of long bursts and their spectral lags (Norris et al. 2000). 
Various authors have studied this relation using arbitrary observer- 
frame energy-bands of various instruments (Ukwatta et al. 2010a; 
Hakkila et al. 2008; Schaefer 2007; Gehrels et al. 2006; Norris 
2002). These investigations support the existence of the relation, 
however with considerable scatter in the extracted results. Re- 
cently, Margutti et al. (2010) investigated spectral lags of X-ray 
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flares and found that X-ray flares of long GRBs also exhibit the 
lag-luminosity correlation observed in the prompt emission. 

The spectral lag is defined as the difference in time of arrival 
of high and low energy photons and is considered to be positive 
when the high-energy photons arrive earlier than the low energy 
ones. Typically the spectral lag is extracted between two arbitrary 
energy bands in the observer frame. However, because of the red- 
shift dependence of GRBs, these two energy bands can correspond 
to a different pair of energy bands in the GRB source-frame, thus 
potentially introducing an arbitrary energy dependence to the ex- 
tracted spectral lag. 

In order to explore whether the lag-luminosity relation is in- 
trinsic to the GRB, it is preferable to extract spectral lags in the 
source frame as opposed to the observer frame. At least two cor- 
rections are needed to accomplish this: 1) Correct for the time dila- 
tion effect (z-correction), and 2) Take into account the fact that for 
GRBs with various redshifts, observed energy bands correspond 
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to different energy bands at the GRB source-frame (K-correction; 
Gehrels et al. (2006)). 

The first correction is straightforward and is achieved by mul- 
tiplying the extracted lag value (in the observer-frame) by (1 + 
z) _1 . The second correction, on the other hand, is not so straight- 
forward. Gehrels et al. (2006) attempted to approximately correct 
the spectral lag by multiplying the lag value (in the observer- 
frame) by (1 + z) 0,33 . We note here that this correction is based 
on the assumption that the spectral lag is proportional to the pulse 
width and that the pulse width itself is proportional to the energy 
(Zhang et al. 2009; Fenimore et al. 1995). These approximations 
depend on clearly identifying corresponding pulses in the light 
curves of each energy band, and may be of limited validity for a 
large fraction of GRBs in which the light curves are dominated by 
overlapping multi-pulse structures. 

Using a sample of 31 Swift GRBs, Ukwatta etal. (2010a) 
(hereafter U10) found that the correlation coefficient improves sig- 
nificantly after the z-correction is applied. Flowever, this correlation 
does not improve further after the application of the K-correction 
as defined by Gehrels et al. (2006). 

An alternative is to make the K-correction by choosing two 
appropriate energy bands fixed in the GRB source-frame and pro- 
jecting these bands into the observer-frame using the relation 
^observer = i?source/(l + z). Ukwatta et al. (2010b) used this 
method for the first time to investigate the lag-luminosity relation 
in the source-frame of the GRB. They selected two source-frame 
energy bands (100 - 200 keV and 300 - 400 keV) and used back- 
ground subtracted as well as non-background subtracted Swift 
data to extract lags. Non-background subtracted data were used to 
improve the signal-to-noise ratio for weak bursts. They found that 
the source-frame relation seems a bit tighter, but with a slope con- 
sistent with previous studies. Arimoto et al. (2010) also looked at a 
limited sample of HETE-II bursts (8 GRBs) both in the observer- 
frame and the source-frame and concluded that there is no signif- 
icant effect from the redshift. However, the redshift distribution of 
their burst sample is very narrow and peaks around one. In contrast 
to Ukwatta et al. (2010b), in this study we used only background 
subtracted data and measured the lag between source- frame energy 
bands 100 - 150 keV and 200 - 250 keV (the reason for selecting 
these particular energy bands is described in § 2) for a sample of 43 
Swift bursts with spectroscopic redshifts. 

In this work we have investigated only long GRBs, i.e., bursts 
with duration greater than ~ 2 seconds. It is rather difficult to 
test the Lag-Liso relation effectively for short GRBs due to a lack 
of spectroscopically measured redshifts. None of the short bursts 
detected so far have any redshift measurements obtained from 
a spectroscopic analysis of their optical afterglow. Moreover, it 
has been shown that short GRBs have either small or negligible 
lags (Norris & Bonnell 2006; Zhang et al. 2006). According to the 
Lag-Li S o relation, these small lag values imply that short bursts to 
be highly luminous. However, based on the redshift measurements 
of their host galaxies we can show that short GRBs are generally 
less luminous than long bursts. Hence short bursts seem to not fol- 
low the lag-luminosity relation (Gehrels et al. (2006)). 

The structure of this paper is the following: In § 2 we dis- 
cuss briefly our methodology for extracting spectral lags. In § 3 we 
present our results for a sample of 43 Swift GRBs. We discuss our 
results with two candidate models in § 4. Finally, in the last section 
(§ 5) we summarize our results and conclusions. Throughout this 
paper, the quoted uncertainties are at the 68% confidence level. 


2 METHODOLOGY 


The Swift Burst Alert Telescope (BAT) is a highly sensitive 
instrument using a coded-mask aperture (Barthelmy et al. 2005). 
BAT uses the shadow pattern resulting from the coded mask to fa- 
cilitate localization of the source. When a gamma-ray source illu- 
minates the coded mask, it casts a shadow onto a position-sensitive 
detector. The shadow cast depends on the position of the gamma- 
ray source on the sky. If one knows the tile pattern in the coded 
mask and the geometry of the detector, it is possible to calculate 
the shadow patterns created by all possible points in the sky using a 
ray-tracing algorithm. Hence, by correlating the observed shadow, 
with the pre-calculated shadow, one can find the location of the 
source. However, each detector can be illuminated by many sources 
and a given source can illuminate many detectors. Hence, in or- 
der to disentangle each sky position, special algorithms have been 
developed and integrated in to the data analysis software by the 
Swift BAT team. 

To generate background-subtracted light curves we used a pro- 
cess called mask weighting. The mask weighting assigns a ray- 
traced shadow value for each individual event, which then en- 
ables the user to calculate light curves or spectra. We used the 
batmaskwtevt and batbinevt tasks in FTOOLS to gener- 
ate mask weighted, background-subtracted light curves, for various 
observer-frame energy bands, as shown in Table 1 . These are the 
energy bands that correspond to fixed energy bands in the source- 
frame i.e. 100 — 150 and 200 — 250 keV. These particular en- 
ergy bands were selected so that after transforming to the observer- 
frame they lie in the detectable energy range of the Swift BAT 
instrument (see Fig. 1). Even though the BAT can detect photons 
up to 350 keV, we limited the upper-boundary to 200 keV in the 
observer-frame. This is because the mask-weighted effective area 
of the detector falls rapidly after 200 keV and as a result the contri- 
bution to the light curve from energies greater than ~200 keV (in 
observer-frame) is negligible (Sakamoto et al. 201 1). 

The spectral lags were extracted using the improved cross- 
correlation function (CCF) analysis method described in U10. In 
this method, the spectral lag is defined as the time delay corre- 
sponding to the global maximum of the cross-correlation function. 
The CCF with a delay index d is defined as, 


CCF(d, x, y) 


E m'm(N,N—d) 

i=max(i,i-d) Xi yi+d 


VEi x l E* yf 


(i) 


where Xi and yi are two sets of time-sequenced data spread over 
N bins. The time delay is obtained by multiplying d by the time 
bin size of the light curves. A Gaussian curve was fitted to the CCF 
(plotted as a function of time delay) to extract the spectral lag. The 
uncertainty in the spectral lag is obtained by simulating 1,000 light 
curves using the Monte Carlo technique (see U10 for more details). 

The isotropic peak luminosity (Li so ) and its uncertainty for 
each GRB is obtained using the method described in U10. In 
essence, a typical GRB spectrum can be described by the Band 
function (Band et al. 1993), for the photon flux per unit photon en- 
ergy using 


N(E) = 



“ e ~( 2 + a ) E / E P k 

) p r («- ffl-Spk J/3-c) , 

) L(2+a)100keVJ C ’ 


E < (i^Pk 


( 2 ) 


which has four model parameters: the amplitude (A), the low- 
energy spectral index (a), the high-energy spectral index (/3) and 
the peak (E p k) of E 2 N(E) spectrum (also called the vF u spec- 
trum, apart from a factor of Planck’s constant). Using these spectral 
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Figure 1. Fixed energy bands at the GRB source-frame are projected to various energy bands at the observer-frame, depending on the redshift. 


parameters the observed peak flux can be calculated for the source- 
frame energy-range E\ = 1.0 keV to E 2 = 10, 000 keV using 

l‘E 2 /(l+z) 

fobs = / N(E)EdE. (3) 

JE 1 /(l+z) 

The isotropic peak luminosity is defined by 

L iso = 47rd| / obs (4) 

where cLl is the luminosity distance: 

d L = f . = (5) 

Jo \/Dm(1 4~ -{- FIl 

For the current universe we take Q.m = 0.27, Q.l = 0.73 and 
the Hubble constant Ho to be 70 (kms _1 )/Mpc (Komatsu et al. 
2009). For more details of the Li so calculation see U10. 


3 RESULTS 

We employed an additional 12 long bursts to the GRB sample (31 
GRBs) that was used in U10, which increased the total sample to 
43. This sample has redshifts ranging from 0.346 (GRB 061021) to 
5.464 (GRB 060927) with an average redshift of ~2.0. The spectral 
information for the additional 12 bursts used in this paper is given in 
Table 2. The calculated peak isotropic luminosities, spanning three 
orders of magnitude, are given in U10 and Table 2. 

By choosing appropriate energy bands in the observer- 
frame (according to the redshift of each burst), we extracted 
mask-weighted background-subtracted light curves for the se- 
lected source-frame energy bands 100-150 and 200-250 keV. The 
observer-frame energy bands used for each burst are shown in Ta- 
ble 1 . Note that the energy gap between the mid-point of the two 
source- frame energy bands is fixed at 100 keV whereas in the 
observer-frame, as expected, this gap varies depending on the red- 
shift of each burst (see the Table 1). For example, in GRB 060927, 
this gap is 16 keV and in GRB 061021, it is 75 keV. This is in con- 
trast to the spectral lag extractions performed in the observer- frame 
where this gap is treated as a constant. 

The extracted spectral lags for the source-frame energy bands 


100-150 and 200-250 keV are listed in Table 3. The Swift BAT 
trigger ID, the segment of the light curve used for the lag extraction 
( T + Xs and T + Xe, T is the trigger time), the time binning of 
the light curve, and the Gaussian curve fitting range of the CCF 
vs time delay plot (with start time, and end time denoted as LS 
and LE respectively) are also given in Table 3. Of 43 bursts in the 
sample there are 24 bursts which have lags greater than zero. The 
remaining 19 bursts either have lags consistent with zero (16 bursts) 
or negative values (3 bursts). 

For the 24 bursts which have positive lags with significance 
l -a or greater (see Table 3), we find that the redshift corrected lag 
is anti-correlated with Li ao . The correlation coefficient for this re- 
lation is -0.82 ± 0.05 with a chance probability of ~ 5.54 x 10~ 5 . 
The extracted correlation coefficient is significantly higher than the 
correlation coefficient (averaged over the six combinations of stan- 
dard BAT energy channels) of ~ —0.68 reported in U10. Various 
correlation coefficients of the relation are shown in Table 4, where 
uncertainties in the correlation coefficients were obtained through 
a Monte Carlo simulation utilizing uncertainties in Li so and the lag 
values. The null probability that the correlation occurs due to ran- 
dom chance is also given for each coefficient type. 

Fig. 2 shows a log- log plot of isotropic peak luminosity vs 
redshift-corrected spectral lag. The solid line shows the following 
best-fit power-law curve: 


iog ( = (54 - 7 ± °- 4) ■ (i - 2 ± °- 2) iog L T+r } - (6) 


Since there is considerable scatter, the uncertainties of the 
fit parameters are multiplied by a factor of ^/x^/ndf = 
\J 84.36/22 — 1.96. The dash lines indicate the estimated 1 -<t 
confidence level, which is obtained from the cumulative fraction of 
the residual distribution taken from 16% to 84%. 

The best-fit power-law index (— 1.2 ± 0.2) is consistent with 
observer- frame results obtained by Norris et al. (2000) (~ —1.14) 
and the average power-law index of — 1.4 ± 0.3 reported in U10. 
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Figure 2. The spectral lags between the source-frame energy range bands 100 — 150 keV and 200 — 250 keV and the isotropic peak luminosity are plotted 
in a log-log plot. 








Figure 3. All combinations of fixed observer-frame energy channel (canonical BAT energy bands: channel 1 (15-25 keV), 2 (25-50 keV), 3 (50-100 keV) 
and 4 (100-200 keV)) spectral lag values as a function of fixed source-frame energy channel (between 100-150 keV and 200-250 keV) lag values. Black and 
red data points and labels corresponds to redshift uncorrected and corrected cases respectively. The blue dashed line corresponds to the equality line of the two 
parameters in each panel. 
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Table 1 . The observer-frame energy bands and energy gaps (the energy difference between the mid-points 
of energy bands) for bursts in the sample. 


GRB 

Redshift 

Low Energy Band (keV) 

High Energy Band (keV) 

Energy Gap (keV) 

GRB050401 

2.899 1 

26-38 

51-64 

26 

GRB050603 

2.821 2 

26-39 

52-65 

26 

GRB050922C 

2.199 3 

31-47 

63-78 

32 

GRB05 1 1 1 1 

1.549 4 * * 

39-59 

78-98 

39 

GRB060206 

4.056 s 

20-30 

40-49 

20 

GRB060210 

3.913® 

20-31 

41-51 

21 

GRB060418 

1.490 7 

40-60 

80-100 

40 

GRB060904B 

0.703 s 

59-88 

117-147 

59 

GRB060908 

1.884 9 

35-52 

69-87 

35 

GRB060927 

5.464 10 

15-23 

31-39 

16 

GRB061007 

1.262 11 

44-66 

88-111 

45 

GRB061021 

0.346 12 

74-111 

149-186 

75 

GRB061 121 

1.315 13 

43-65 

86-108 

43 

GRB070306 

1.496 14 

40-60 

80-100 

40 

GRB071010B 

0.947 15 

51-77 

103-128 

52 

GRB071020 

2.145 16 

32-48 

64-79 

32 

GRB080319B 

0.937 17 

52-77 

103-129 

52 

GRB080319C 

1.949 18 

34-51 

68-85 

34 

GRB08041 1 

1.030 19 

49-74 

99-123 

50 

GRB080413A 

2.433 20 

29-44 

58-73 

29 

GRB080413B 

1.101 21 

48-71 

95-119 

48 

GRB080430 

0.767 22 

57-85 

113-141 

56 

GRB080603B 

2.689 23 

27-41 

54-68 

27 

GRB080605 

1.640 24 

38-57 

76-95 

38 

GRB080607 

3.036 25 

25-37 

50-62 

25 

GRB080721 

2.591 26 

28-42 

56-70 

28 

GRB080916A 

0.689 27 

59-89 

118-148 

59 

GRB081222 

2.770 28 

27-40 

53-66 

26 

GRB090424 

0.544 29 

65-97 

130-162 

65 

GRB090618 

0.540 30 

65-97 

130-162 

65 

GRB090715B 

3.000 31 

25-38 

50-63 

25 

GRB090812 

2.452 32 

29-43 

58-72 

29 

GRB090926B 

1.240 33 

45-67 

89-112 

45 

GRB091018 

0.971 34 

51-76 

101-127 

51 

GRB091020 

1.710 35 

37-55 

74-92 

37 

GRB091024 

1.091 36 

48-72 

96-120 

48 

GRB091029 

2.752 37 

27-40 

53-67 

27 

GRB091208B 

1.063 38 

48-73 

97-121 

49 

GRB 100621 A 

0.542 39 

65-97 

130-162 

65 

GRB100814A 

1.440 40 

41-61 

82-102 

41 

GRB 10081 6A 

0.800 41 

56-83 

111-139 

56 

GRB 1 00906A 

1.727 42 

37-55 

73-92 

37 

GRB1 10205A 

2.220 43 

31-47 

62-78 

31 

GRB1 10213A 

1.460 44 

41-61 

81-102 

41 


(1) Watson et al. (2006); (2) Berger & Becker (2005b); (3) Piranomonte et al. (2008); (4) Penprase et al. 
(2006); (5) Fynbo et al. (2009); (6) Fynbo et al. (2009); (7) Prochaska et al. (2006); (8) Fynbo et al. (2009); 
(9) Fynbo et al. (2009); (10) Fynbo et al. (2009); (11) Fynbo et al. (2009); (12) Fynbo et al. (2009); (13) 
Fynbo et al. (2009); (14) Jaunsen et al. (2008); (15) Cenko et al. (2007); (16) Jakobsson et al. (2007); (17) 
D’Elia et al. (2009); (18) Fynbo et al. (2009); (19) Fynbo etal. (2009); (20) Fynbo et al. (2009); (21) 
Fynbo et al. (2009); (22) Cucchiara & Fox (2008); (23) Fynbo et al. (2009); (24) Fynbo et al. (2009); (25) 
Prochaska et al. (2009); (26) Fynbo et al. (2009); (27) Fynbo et al. (2009); (28) Cucchiara et al. (2008); (29) 
Chomock et al. (2009); (30) Cenko et al. (2009); (31) Wiersema et al. (2009); (32) de Ugarte Postigo et al. 
(2009); (33) Fynbo et al. (2009); (34) Chen et al. (2009); (35) Xu et al. (2009); (36) Cucchiara et al. 
(2009); (37) Chomock et al. (2009); (38) Wiersema et al. (2009); (39) Milvang-Jensen et al. (2010); (40) 
O’Meara et al. (2010); (41) Tanvir et al. (2010); (42) Cenko et al. (2011); (43) Milne & Cenko (2011). 


4 DISCUSSION 

4.1 Spectral Lags: Observer-frame versus Source-frame 

U10 extracted spectral lags in fixed energy bands in the observer- 

frame and in this work for the same sample of 31 bursts we ex- 


tracted lags in fixed energy bands in the source-frame. In the 
observer-frame case, there are four energy channels (canonical 
BAT energy bands: channel 1 (15-25 keV), 2 (25-50 keV), 3 (50- 
100 keV) and 4 (100-200 keV)), thus six lag extractions per burst. 
It is interesting to study to what degree these different lags cor- 
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Table 2. GRB redshift and spectral information. Note that uncertainties of parameters that are reported with 90% confidence level have been reduced to 
l-cr level for consistency. 


GRB 

Peak Flux a 

E pk b 

a 

/3 C 

L iso erg/s 


Reference 

GRB090812 

3.60 ±0.13 

C79 + 99 
°' z -156 

— 1 03" 1 ” 0 ' 04 
-L.uo_o 04 

—2 ^n+O- 16 
z.ou _ 0 16 

(7-86 tifr) 

X 

10 s2 

PaTShin et al. (2009); Baumgartner et al. (2009) 

GRB090926B 

3.20 ±0.19 

91+1 

—0 13+ 0 - 04 
u ‘ iO -0.04 

—2 36+0- 31 
z.oo_o 31 

(5.22+g;“) 

X 

10 51 

Briggs (2009); Baumgartner et al. (2009) 

GRB09 1018 

10.30 ±0.25 

28+g° 

i ro+0.24 
-L.oo_o 37 

—2 44 + 0,15 
0.15 

(6.96±J;™) 

X 

10 51 

Golenetskii et al. (2009); Markwardt et al. (2009) 

GRB09 1 020 

4.20 ±0.19 

+ 1 

b- 

—0 20” 1 ” 0 ' 25 
U.zu_o 25 

—1 70+00 1 

(2.8l±g;ig) 

X 

10 52 

Chaplin (2009); Palmer et al. (2009) 

GRB091024 

2.00 ±0.19 

03 

O 

0 

1 + 

O O 
O O 

— 1 lo - * -013 

A * 1U — 0.13 

—2 36+0- 31 
z,c>o -0.31 

(5.56l 2 ;«) 

X 

10 51 

Golenetskii et al. (2009); Sakamoto et al. (2009) 

GRB091029 

1.80 ±0.06 

6 i + 10 

D± -10 

— 1 46" 1 ” 0 ' 17 
0.17 

—2 36+ 0,31 
^.ou_o 31 

(X.6T±g;“) 

X 

10 52 

Barthelmy et al. (2009) 

GRB091208B 

15.20 ±0.63 

124 + 2 

_i 44+0-04 
A *”— 0.04 

r\ oq + 0.29 
0.12 

(l- 68 ±g; 6 |) 

X 

10 52 

McBreen (2009); Baumgartner et al. (2009) 

GRB 100621 A 

12.80 ±0.19 

95l?r 

—1 7 o+ 0,08 
A * /u _0.08 

-2 45 + 1 - 44 

(2-55t°j®) 

X 

10 51 

Golenetskii et al. (2010); Ukwatta et al. (2010) 

GRB100814A 

2.50 ±0.13 

106+g 

—0 64 +0 ' 08 
u-o4:_o 09 

9 r» 9 + 0.08 
^• u ^-0.06 

(8-27 ±S:Jg) 

X 

10 s1 

von Kienlin (2010); Krimm et al. (2010) 

GRB100906A 

10.10 ±0.25 

180+28 

—1 io +006 
±1U - 0.06 

—2 20”*” 0 ' 19 
^.^u_o 13 

(4-90l+ 2 3) 

X 

10 52 

Golenetskii et al. (2010); Barthelmy et al. (2010) 

GRB 11 0205 A 

3.60 ±0.13 

999+46 

ZZZ -46 

1 ^9+0-09 
A '^ z — 0.09 

—2 36+ 0,31 

Z.OO-o 31 

(2-78 ±g;®5) 

X 

10 52 

Golenetskii et al. (2011) 

GRB1 10213A 

1.60 ±0.38 

98^5 

— 1 44+° 03 
A,f ”— 0.03 

—2 36+0 31 

^.ou _0 31 

( 3 . 53 +^ 53 ) 

X 

10 51 

Foley (2011); Barthelmy et al. (2011) 


“ 1-second peak photon flux measured in photons cm -2 s — 1 in the energy range 15 — 150 keV. 
b Peak energy, E pk , is given in keV. 

c Values in brackets indicates estimated high-energy photon index, /3, which is the mean value of the BATSE /3 distribution (Kaneko et al. 2006; 
Sakamoto et al. 2009). 


Table 4. Correlation coefficients of the lag-luminosity relation. 


Coefficient Type 

Correlation Coefficient 

Null Probability 

Pearson’s r 

- 0.82±0.05 

5.54 X 10 -5 

Spearman’s r s 

- 0.70±0.06 

1.49 x 10 -4 

Kendall’s r 

- 0.50±0.05 

6.63 x 10 -4 


relate with source-frame lags (between fixed source-frame energy 
channels 100-150 keV and 200-250 keV). In Fig. 3 we show all 
combinations of observer-frame lags as a function of source-frame 
lags. The red data points show lags with the time-dilation correction 
due to cosmological redshift and black data points show lags with- 
out the time-dilation correction. From Fig. 3 it is clear that all plots 
show some correlation both in the time-dilation corrected (shown 
in red) and time-dilation uncorrected (shown in black) cases. We 
note that the correlation coefficients are greater than 0.5 in time- 
dilation uncorrected cases where BAT channel 1 is involved in the 
lag extraction. In the time-dilation corrected case all plots show cor- 
relation coefficients greater than 0.5 except for the Lag43 plot. De- 
spite these moderate correlation coefficients, the large scatter seen 
in these plots indicate that the observer-frame lag does not directly 
represent the source-frame lag. 


4.2 Fag- /.is,, Relation: Observer-frame versus Source-frame 

There are two important changes in the Lag-Luminosity relation 
which may occur when going from fixed observer-frame energy 
bands to fixed source-frame energy bands: A change in the power- 
law index, and a change in the dispersion of the data measured by 
the correlation coefficient. Table 5 summarizes these two parame- 
ters for various energy bands both in the observer-frame and in the 
source- frame. 

In the observer-frame the power-law index varies from ~ 0.6 
to ~ 1.8, with mean around 1.3. In the source-frame the index 
changes from 0.9 to 1 .23 with a mean of ~ 1 . 1 . Meanwhile, the cor- 
relation coefficient varies from 0.60 to 0.79 in the observer-frame 
and in the source-frame it changes from 0.76 to 0.90. Hence, ac- 


cording to Table 5, the source-frame Lag-Li so relation seems to be 
tighter than the observer- frame case with a slope closer to one. 


4.3 Spectral Lag - E p k Relation 

Now we investigate the relation between source-frame spectral lag 
and source-frame average peak energy (i7 p k(l + 2 )) of the burst 
spectrum. In Fig. 4, we plotted E p k(l + z) as a function of source- 
frame lags. There is a correlation between these two parameters 
with a correlation coefficient of —0.57 ± 0.14. Various correlation 
coefficients of the relation are shown in Table 6, with uncertainties 
and null probabilities. 

The best-fit is shown as a dashed line in Fig. 4, yielding the 
following relation between .E p k(l + z ) and Lag/( 1 + z): 

log - (3.7±0.1)— (0.56±0.06) log .(,) 

The uncertainties in the fitted parameters are expressed with 
the factor of \Jx 2 / n df = +30.71/22 « 1.18. 

According to equation (6), Li so oc ( Lag/{1 + z ))~ 1,2 . 
From the Yonetoku relation we know that Li so oc (_E p k(l ± 
z)) 2 ' 0 (Yonetoku et al. 2004). Hence, from these two relations we 
expect to see a correlation between i? p k(l + z) and Lag /(I + z) 
such as i7 p k(l + z) oc (Lag/( 1 ± z)) -0 ’ 6 . 

The best fit slope of 0.56±0.06 is consistent with the expected 
slope of ~ 0.6 based on the source-frame lag-luminosity and the 
Yonetoku relation. However, note that the correlation coefficient is 
significantly smaller than the coefficient for the lag-luminosity rela- 
tion. This lower degree of correlation may be suggestive of bright- 
ness and detector related selection effects that have been noted in 
the literature (Butler et al. 2007) for the Yonetoku relation. 


4.4 Some Models for Spectral Lags 

U 1 0 and this work have provided more evidence for the existence of 
the lag-luminosity relation based on a sample of Swift BAT GRBs 
with measured spectroscopic redshifts. This analysis calls for a 
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Table 3. Source-frame spectral lag values of long duration Swift BAT GRBs 


GRB 

Trigger ID 

T + X s (s) 

T + X e (s) 

Bin Size (ms) 

LS (s) 

LE (s) 

Lag Value (ms) 

Significance 

GRB05040 1 

113120 

23.03 

29.43 

64 

-2.00 

2.00 

310±145 

2.14 

GRB050603 

131560 

-3.83 

3.08 

16 

-0.40 

0.40 

-16±21 

-0.76 

GRB050922C 

156467 

-2.70 

2.94 

16 

-1.00 

1.00 

136±68 

2.00 

GRB05 1111 

163438 

-6.96 

28.62 

64 

-4.00 

4.00 

333±251 

1.33 

GRB060206 

180455 

-1.29 

8.18 

16 

-2.00 

2.00 

86± 111 

0.77 

GRB060210 

180977 

-3.37 

5.08 

128 

-4.00 

4.00 

658±259 

2.54 

GRB0604 1 8 

205851 

-7.66 

33.04 

64 

-2.00 

2.00 

-110±106 

-1.04 

GRB060904B 

228006 

-1.97 

10.32 

512 

-6.00 

6.00 

124±436 

0.28 

GRB060908 

228581 

-10.91 

3.68 

32 

-2.00 

2.00 

783=124 

0.63 

GRB060927 

231362 

-1.69 

8.04 

32 

-1.00 

1.00 

183=75 

0.24 

GRB061007 

232683 

23.86 

65.08 

4 

-0.20 

0.20 

523=22 

2.36 

GRB06 1 02 1 

234905 

-0.46 

14.64 

512 

-4.00 

4.00 

-4303=975 

-0.44 

GRB061 121 

239899 

60.44 

80.66 

4 

-0.20 

0.20 

22±10 

2.20 

GRB070306 

263361 

90.00 

118.42 

32 

-4.00 

2.00 

-362±247 

-1.47 

GRB071010B 

293795 

-1.70 

17.24 

64 

-2.00 

2.00 

4043=159 

2.54 

GRB071020 

294835 

-3.22 

1.14 

4 

-0.20 

0.40 

35±13 

2.69 

GRB080319B 

306757 

-2.85 

57.57 

4 

-0.10 

0.14 

23±6 

3.83 

GRB080319C 

306778 

-0.77 

13.31 

32 

-1.00 

1.00 

174±91 

1.91 

GRB0804 1 1 

309010 

38.46 

48.45 

4 

-0.50 

0.50 

116±25 

4.64 

GRB080413A 

309096 

-0.42 

9.05 

8 

-1.00 

1.00 

1073=59 

1.81 

GRB080413B 

309111 

-1.44 

4.96 

32 

-1.00 

1.00 

1153=50 

2.30 

GRB080430 

310613 

-1.24 

12.84 

256 

-4.00 

4.00 

9 1 ±43 1 

0.21 

GRB080603B 

313087 

-0.54 

5.10 

16 

-1.00 

1.00 

5±59 

0.08 

GRB080605 

313299 

-5.46 

15.53 

8 

-0.20 

0.20 

35±18 

1.94 

GRB080607 

313417 

-6.13 

12.05 

8 

-0.50 

0.50 

26±30 

0.87 

GRB080721 

317508 

-3.39 

8.64 

64 

-2.00 

2.00 

-863=110 

-0.78 

GRB080916A 

324895 

-2.66 

39.58 

128 

-2.00 

4.00 

5853=214 

2.73 

GRB081222 

337914 

-0.80 

15.58 

4 

-1.00 

1.00 

2273=51 

4.45 

GRB090424 

350311 

-0.94 

4.95 

16 

-0.20 

0.20 

14±14 

1.00 

GRB0906 1 8 

355083 

46.01 

135.35 

8 

-2.00 

2.00 

267±72 

3.71 

GRB090715B 

357512 

-4.80 

21.06 

16 

-2.00 

3.00 

2753=155 

1.77 

GRB090812 

359711 

-6.93 

41.20 

256 

-6.00 

6.00 

-223=202 

-0.11 

GRB090926B 

370791 

-22.00 

36.00 

512 

-10.00 

8.00 

746±627 

1.19 

GRB091018 

373172 

-0.28 

2.92 

64 

-2.00 

1.00 

1433=297 

0.48 

GRB091020 

373458 

-2.54 

13.84 

128 

-3.00 

2.00 

-187±177 

-1.06 

GRB091024 

373674 

-9.58 

27.29 

512 

-10.00 

10.00 

9123=604 

1.51 

GRB09 1 029 

374210 

-4.03 

38.98 

256 

-10.00 

10.00 

-1 12±395 

-0.28 

GRB091208B 

378559 

7.66 

10.61 

64 

-1.00 

1.00 

105±66 

1.59 

GRB 100621 A 

425151 

-6.79 

40.31 

256 

-3.00 

3.00 

1199±311 

3.86 

GRB100814A 

431605 

-4.40 

29.39 

256 

-4.00 

4.00 

862±147 

5.86 

GRB100906A 

433509 

-1.49 

26.16 

128 

-2.00 

2.00 

1053=79 

1.33 

GRB 11 0205 A 

444643 

118.89 

293.99 

64 

-1.00 

1.00 

-29±52 

-0.56 

GRB1 10213A 

445414 

-3.42 

5.29 

512 

-3.00 

3.50 

6023=746 

0.81 


Table 6. Correlation coefficients of the lag-E p ^ relation. 


Coefficient Type 

Correlation Coefficient 

Null Probability 

Pearson’s r 

- 0.573=0. 14 

4.83 X 10“ 3 

Spearman’s r s 

- 0.503=0. 12 

1.36 x 10" 2 

Kendall’s r 

- 0.37±0.14 

1.18 x 10 -2 


physical interpretation for spectral lag and a lag- luminosity rela- 
tion. In the literature, several possible interpretations have been dis- 
cussed (Dermer 1998; Salmonson 2000; Ioka & Nakamura 2001; 
Kocevski & Liang 2003; Schaefer 2004; Qin et al. 2004; Ryde 
2005; Shen et al. 2005; Lu et al. 2006; Peng et al. 2011). 

One proposed explanation for the observed spectral lag is the 
spectral evolution during the prompt phase of the GRB (Dermer 
1998; Kocevski & Liang 2003; Ryde 2005). Due to cooling effects, 


.Epk moves to a lower energy channel after some characteristic 
time. When the peak energy (E p k) moves from a higher energy 
band to a lower energy band, the temporal peak of the light curve 
also moves from a higher energy band to a lower one, which results 
in the observed spectral lag. In a recent study, Peng et al. (2011) 
suggest that spectral evolution can be invoked to explain both pos- 
itive and negative spectral lags. Hard-to-soft evolution of the spec- 
trum produces positive spectral lags while soft-to-hard evolution 
would lead to negative lags. In addition, these authors also suggest 
that soft-to-hard-to-soft evolution may produce negative lags. 

A schematic diagram showing a hard-to-soft scenario is de- 
picted in Fig. 5. Initially, _E p k of the spectrum is in the high-energy 
band, which results in a pulse in the light curve of the high en- 
ergy band. Then E p k moves to the lower energy band resulting in 
a pulse in the low-energy light curve. The temporal difference be- 
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Table 5. Observer-frame and source-frame slopes and correlation coefficients of the lag-I/i so relation. Conservative 10% 
uncertainty is assumed for cases without uncertainties. 


Energy Bands 

Frame 

Slope 

Correlation Coefficient 

Number of GRBs 

Reference 

(0.3-1), (3-10) keV 

Observer 

0.95±0.23 

- 

9 

Margutti et al. (2010) 

(6-25), (50-400) keV 

Observer 

1.163=0.07 

- 0 79+ 0 - 16 
U ‘ ''— 0.05 

8 

Arimoto et al. (2010) 

(15-25), (25-50) keV 

Observer 

1.43=0.1 

- 0.633=0.06 

21 

U10 

(15-25), (50-100) keV 

Observer 

1.53=0.1 

- 0.603=0.06 

28 

U10 

(15-25), (100-200) keV 

Observer 

1.83=0.1 

- 0.673=0.07 

27 

U10 

(25-50), (50-100) keV 

Observer 

1.23=0.1 

- 0.663=0.07 

27 

U10 

(25-50), (100-200) keV 

Observer 

1.43=0.1 

- 0.753=0.07 

25 

U10 

(25-50), (100-300) keV 

Observer 

1.143=0.1 

- 

6 

Norris et al. (2000) 

(25-50), (100-300) keV 

Observer 

0.623=0.04 

-0.723=0.07 

6 

Hakkila et al. (2008) 

(50-100), (100-200) keV 

Observer 

1.43=0.1 

- 0.773=0.08 

22 

U10 

(20-100), (100-500) keV 

Source 

1.233=0.07 

- 0 90 +0 ' 12 
u y - 0.02 

8 

Arimoto et al. (2010) 

(100-200), (300-400) keV 

Source 

0.93=0.1 

- 0.763=0.06 

22 

Ukwatta et al. (2010b) 

(100-150), (200-250) keV 

Source 

1.23=0.2 

- 0.823=0.05 

24 

This work 



Lag/(1+z) msec 

Figure 4. The source-frame peak energy (E p k(l + 2 )) versus source-frame spectral lags. The energy bands, 100 — 150 keV and 200 — 250 keV, corresponding 
to the lag extractions are shown in hashed red bands on the plot. 


tween the two pulses in the light curves would then be a measure 
of the cooling time scale of the spectrum. 

If this were the only process that caused the lag then in a sim- 
ple picture one would expect the source-frame average E p y to lie 
within the two energy bands in question. According to Fig. 4, for 
the majority of bursts the source-frame _E p k lies outside the energy 
band 100 — 250 keV, indicating that the simple spectral evolution 
scenario described above may not be the dominant process respon- 
sible for the observed lags. However, it is worth noting that a pulse 
in a specific energy band may not always mean that the E p k is also 
within that energy band. There are other issues associated with this 
model: 1) the calculated cooling times based on simple synchrotron 
models are, in general, relatively small compared to the extracted 
lags, and 2) short bursts which exhibit considerable spectral evolu- 
tion do not show significant lags. 


Another model that purports to explain spectral lags is based 
on the curvature effect, i.e., a kinematics effect due to the observer 
looking at increasingly off-axis annulus areas relative to the line- 
of-sight (Salmonson 2000; Ioka & Nakamura 2001; Dermer 2004; 
Shenetal. 2005; Luetal. 2006). Fig. 6 illustrates how the spec- 
tral lag could arise due to the curvature effect of the shocked shell. 
Due to a smaller Doppler factor and a path difference, the radia- 
tion from shell areas which are further off axis will be softer and 
therefore lead to a lag. As with spectral evolution models, there are 
difficulties associated with the curvature models too. These kine- 
matic models generally predict only positive lags. As can be seen 
from Table 3 some of the measured lags are negative, and therefore 
these lags present a real challenge for the simple curvature models. 

It is possible that spectral lags are caused by multiple mech- 
anisms. Peng et al. (2011) investigated spectral lags caused by in- 
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Figure 5. The time evolution of the E p k across energy bands may cause the observed spectral lags in GRBs. 


Hard 



Figure 6. Spectral lags could arise due to the curvature effect of the shocked shell. At the source, the relativistically expanding shell emits identical pulses from 
all latitudes. However, when the photons reach the detector, on-axis photons get boosted to higher energy (hard). Meanwhile off-axis photons get relatively 
smaller boost and travel longer to reach the detector. Thus these photons are softer and arrive later than the on-axis photons. 


trinsic spectral evolution and the curvature effect combined. They 
showed that the curvature effect always tends to increase the ob- 
served spectral lag in the positive direction. Even for cases with 
soft-to-hard spectral evolution, when the curvature effect is intro- 
duced lags become positive. Hence they predict that the majority of 
measured spectral lags should be positive, which is consistent with 
the findings of this work and U10. 


5 SUMMARY AND CONCLUSION 


spectral lag (r) and isotropic luminosity over those previously re- 
ported with the following relation: 


t/ ( ms) ^ L iso /(ergs *) 
1 + 2 “ 10 54 - 7 


( 8 ) 


We also find a modest correlation between the source-frame 
spectral lag and the peak energy of the burst, which is given by the 
relation, 


We have investigated the spectral lag between 100 — 150 keV and 
200—250 keV energy bands at the GRB source- frame by projecting 
these bands to the observer- frame. This is a step forward in the 
investigation of lag-luminosity relations since most of the previous 
investigations used arbitrary observer-frame energy bands. 

Our analysis has produced an improved correlation between 


t/ ( ms) ^ Jg p k(l + z)/(keV) 
1 + 2 “ 10 3 - 7 


( 9 ) 


Finally, we mentioned two simple models and noted their lim- 
itations in explaining the observed spectral lags. 
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